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Abstract 

ARGO-YBJ preliminary results of the measurements of the all-particle and light-component (i.e. 
protons and helium) energy spectra between approximately 5 TeV and 5 PeV are reported and 
discussed. The resulting all-particle spectrum (measured in the energy range 100 TeV - 3 PeV) is 
in good agreement with both theoretical parametrizations and previous measurements. The light- 
component (i.e. p + He) has been measured with high resolution up to about 5 PeV. The preliminary 
result is in agreement with direct measurements and then show a clear indication of a bending below 
1 PeV. Improvements of event selection/reconstruction with the full statistics and a complete analysis 
of systematic uncertainties is currently under way. 


1 Introduction 

There is a general consensus that galactic cosmic rays (hereafter CRs) up to the “knee” of the all¬ 
particle spectrum (~ 4 PeV) mainly originate in Supernova Remnants (SNRs). Recent measurements 
carried out by the balloon-borne CREAM experiment mm show that the proton and helium spectra 
from 2.5 to 250 TeV are harder compared to lower energy measurements. As pointed out by several 
authors, the evolution of the proton and helium spectra and their subtle differences could be indications 
of the presence of different populations of CR sources contributing to the overall flux and operating in 
environments with different chemical compositions ei a- Diffusion effects during CR propagation in 
the Galaxy might also play an important role. 

In the knee region (and above) the measurements of the CR primary spectrum are carried out by 
EAS arrays. In this case mass composition studies are extremely difficult and often affected by large 
systematic uncertainties. In the standard picture the average composition at the knee is dominated by 
light elements, and the knee itself is interpreted as the steepening of the p and He spectra [5j. However, 
several experimental results suggest an heavier composition at knee energies ei m 0 is uni mi 

A measurement of the CR primary energy spectrum (all-particle and light-component) in the energy 
range few TeV - 10 PeV is under way with the ARGO-YBJ experiment (for a description of the detector 
and a report of the latest physics results see m)- In order to cover this wide energy range, different 
approaches have been followed: 

- ’digital readout’, based on the RPC readout strip multiplicity, in the 5 TeV - 200 TeV range E3; 

- ’analog readout’, based on the particle density near the shower core, in the 100 TeV-10 PeV range; 

- ’hybrid measurement-’, carried out by ARGO-YBJ and a wide field of view Cherenkov telescope, in 
the 100 TeV - PeV region [H] , 

Preliminary results concerning the all-particle and the light-component (i.e. p+He) spectra 
obtained with the analog readout are summarized in the following. More details on the results obtained 
with the ’hybrid measurement’ are also given in ns ng. 


2 Measurement of the all-particle spectrum 

The measurement of the CR energy spectrum up to 10 PeV is under way exploiting the RPC charge 
readout of the ARGO-YBJ detector which allows studying the structure of the particle density 
distribution in the shower core region up to particle densities of about 10 4 /m 2 Eg. The study 
of the charged particle lateral density function (LDF) at ground is expected to provide information on 
the longitudinal profile of the showers in the atmosphere, that is to estimate their development stage, 
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Figure 1: Primary energy as a function of the 
observed truncated size N p s for simulated showers 
due to different primary nuclei. 


Figure 2: The age parameter s' resulting from the 
fit of the average LDF for simulated showers vs 
the corresponding X max average values. 



Figure 3: All-particle energy spectrum of primary CRs measured by ARGO-YBJ (see text). The 
parametrizations provided by [191 and |5j are shown for comparison. 
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or the so-called age, which is related to X max , the atmospheric depth at which the cascade reaches 
its maximum size. This implies the possibility of selecting showers within given intervals of X max or, 
equivalently, of Xdm, the grammage between the depth of the shower maximum and the detection level. 
For this reason, the combined use of the shower energy and age estimations can ensure a sensitivity 
to the primary mass, thus giving the possibility of selecting a light (p+He) event sample with high 
efficiency. 

Various observables were considered and analyzed in order to find a suitable estimator of the 
primary CR energy. Among them, according to MC simulations, N P 8, the number of particles detected 
within a distance of 8 m from the shower axis, resulted well correlated with energy, not biased by 
the finite detector size and not much affected by shower to shower fluctuations m- Nevertheless, 
as shown in Fig. [I] this truncated size is a mass-dependent energy estimator parameter. In order 
to have a mass-independent parameter we fitted the LDFs of individual showers (up to 10 m from 
the core) event-by-event, for different N p s intervals and different shower initiating primaries, with a 
suitable function to get the shape parameter s' (see DUCES] for details). From these studies we find 
that, for a given primary, the s' value decreases when N p s (i.e. the energy) increases, this being 
due to the observation of younger (deeper) showers at larger energies. Moreover, for a given range 
of N P 8, s' increases going from proton to iron, as a consequence of older (shallower) showers. Both 
dependencies are in agreement with the expectations, the slope s' being correlated with the shower 
age, thus reflecting its development stage. This outcome has two important implications, since the 
measurements of s' and N p s can both (i) help constraining the shower age and (ii) give information on 
the primary particle nature. 

Concerning the first point, we show in Fig. [5] the s' values as obtained from the fit of the average 
LDFs, for each simulated primary type and N p s interval, as a function of the corresponding A max 
average value. As can be seen, the shape parameter s' depends only on the development stage of 
the shower, independently from the nature of the primary particle. That plot expresses an important 
universality of the LDF of detected EAS particles in terms of the lateral shower age. The LDF slope 
s' is a mass-independent estimator of the average Xm ax . Obviously shower-to-shower fluctuations 
introduce unavoidable systematics, whose effects can be anyway quantified and taken into account. 
Another implication is that s' from the LDF fit close to the shower axis, together with the measurement 
of the truncated size N p s, can give information on the primary particle nature, thus making possible 
the study of mass composition and the selection of a light-component data sample. 

Assuming an exponential absorption after the shower maximum, we get a variable linearly 

correlated to the size at the shower maximum, by using N p s and s' measurements for each event and 
simply correcting with an exponential attenuation: N^ ax « N p g ■ exp[(hosec6 — X max (s')) / \ a bs\- A 
suitable choice of the absorption lenght \ a bs (=120 g/cm 1 2 * ) allows to get a parameter correlated 

with primary energy in an almost linear and mass independent way, providing an energy estimator 
with a Log(E/TeV) resolution of 0.10-0.15 (getting better with energy). 

As described in [16], the RPC charge readout system has 8 different and overlapping gain scale 
settings (G0,....,G7 from smaller to larger gains) in order to explore the particle density range «20 
10 4 particles/nr 2 . In this preliminary analysis the results obtained with two gain scales (so-called 
G1 and G4) are presented. Selecting quasi-vertical events (9 < 15°) in terms of the truncated size 
N p g with the described procedure we reconstructed the CR all-particle energy spectrum shown in the 
Fig[3]in the energy range 100 3000 TeV. In the plot a ±14% systematic uncertainty, due to hadronic 

interaction models, selection criteria, unfolding algorithms, and aperture calculation, is shown by the 
shaded area. The statistical uncertainty is shown by the error bars. As can be seen from the figure, 
the two gain scales overlap making us confident about the event selection and the analysis procedure. 
The ARGO-YBJ all-particle spectrum is in fair agreement with the parametrizations provided by m 
and [3], showing evidence of a spectral index change at an energy consistent with the position of the 
knee. As shown in Fig[6] this result is also consistent with previous measurements made by both direct 
and indirect experiments. This is also an important check on the absolute energy scale set for this 
analysis, whose systematic uncertainty has been anyhow conservatively estimated at the level of 10%. 


3 Measurement of the light-component energy spectrum 

The CR light-component energy spectrum has been measured by ARGO-YBJ from about 5 TeV to 
700 TeV [13] [141112) . The energy range is now being extended up to the few PeV region by using the 
RPC charge readout information and three different approaches. 

(1) A selection of events in the s' - N p s space allowing to get a light-component sample of showers 
with a contamination of heavier nuclei less than about 15% (see Fig. [4]). 

(2) A Bayesian unfolding technique similar to that applied to measure the light-component spectrum 

up to 200 TeV [32 132 ED] • 
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Figure 4: Relation between the LDF shape parameter s' and the truncated size N p $ for different nuclei. 
The p+He selection cut is shown by the lines. 


(3) The ARGO-YBJ/WFCTA hybrid measurement [TJ with a different selection procedure which 
increases the aperture of a factor 2.4, thus allowing the extension to larger energies (see jl5j for 
a detailed description of the method and a discussion of the results). 

Preliminary results of the three analyses (and the previous ARGO-YBJ measurement below 
200 TeV) are summarized in Fig[5] The systematic uncertainty on the flux is shown by the shaded area 
and the statistical one by the error bars. A systematic uncertainty on the energy scale at the level of 
5-10% (depending on the analysis) has also been conservatively estimated (not shown in the plots). 
All three different analyses are consistent with low energy (direct) measurements and show a clear 
evidence for a bending at larger energies but below 1 PeV. With respect to a single power-law with 
a spectral index -2.62 the deviation is observed at a level of about 6 s.d. The results obtained with 
the two analysis of RPC charge readout data (label 1 and 2 in the previous list) are in fair agreement. 
They also agree with the ARGO-YBJ/WFCTA hybrid measurement within systematic uncertainties 
and the possible difference in the energy scale. For comparison, the parametrizations of the light- 
component provided by [19] and [3] are shown by the blue and red dashed lines, respectively. A 
Horandel-like spectrum with a modified knee at Zxl PeV (a factor four lower in energy than in the 
original formulation) is also shown for comparison. Finally ARGO-YBJ results are compared to a 
compilation of several other measurements in Fig[6] 


4 Conclusions 

The CR spectrum has been studied by the ARGO-YBJ experiment in a wide energy range (TeVs —¥ 
PeVs) . This study is particularly interesting because not only it allows a better understanding of 
the so called ’knee’ of the energy spectrum and of its origin, but also provides a powerful cross-check 
among very different experimental techniques. The all-particle spectrum (measured in the energy 
range 100 TeV - 3 PeV) is in good agreement with both theoretical parametrizations and previous 
measurements, making us confident about the selection and reconstruction procedures. The light- 
component (i.e. p + He) has been measured with high resolution up to about 5 PeV. The ARGO-YBJ 
preliminary result is in agreement with direct measurements and then show a clear indication of a 
bending below 1 PeV. Improvements of event selection with the full statistics and a complete analysis 
of systematic uncertainties is currently under way. 
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Figure 5: Light (i.e. p+He) component energy spectrum of primary CRs measured by ARGO-YBJ with 
four different analyses (see text). 
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Figure 6: All particle and light (p+He) component energy spectra of primary CR measured by ARGO-YBJ 
and compared to different experimental results. 
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